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Abstract 



A comparison of the rates for ip{2S) — > vr+vr J/tp, J/ip l^l and J/^lJ 
anything is used to determine the J/'0 leptonic branching fractions. The 
results are BiJ/ip e+e") = (5.90 ± 0.05 ± 0.10)% and BiJ/ip n+n') = 
(5.84 lb 0.06 lb 0.10)%, where the first error is statistical and the second is 
systematic. Assuming lepton universality, the leptonic branching fraction of 
the J/ip is B{J/i; £+£-) = (5.87 ±0.04 ±0.09)% per species. This result is 
used to estimate the QCD scale factor A^^^ and the strong coupling constant 

iVi O 

as- 



I. INTRODUCTION 

The branching fractions for the leptonic decays J/if) e+e" {Be) and fj^~^fi~ {B^) are 
basic parameters of the J /if) resonance. They can be used to determine the strong couphng 
constant as or, equivalently, the fundamental scale parameter of QCD, ^jjg [0]. The ratio 
B(./ B^ provides a test of lepton universality. In addition, these branching fractions are used 
to determine the total number of J/ip events in a wide variety of measurements that take 
advantage of the clean experimental J/ip (£ = e or /i) signature. 

The first reported measurement of the leptonic J/ip branching fractions has a precision 
of about 15% and is based on an energy scan across the resonance performed by the Mark-I 
group [0. A subsequent measurement by the Mark-III group is based on a comparison 
of the rates for ip{2S) n^n^J/ip, J/ip ^ and J/ip anything and is, thus, 

independent of the luminosity determination. The Mark-III measurement has a precision of 
4% and is about one (Mark I) standard deviation below the Mark-I result. More recently, 
BES performed an energy scan measurement and obtained results in good agreement with 
Mark-III, but with larger errors [Q. In this paper, we report the results of a measurement 
of the leptonic branching fractions using a sample of ip{2S) decays measured in the BES 
detector at the BEPC storage ring. We apply a technique similar to that used by the 
Mark-III group to a larger data sample. 



II. THE BES DETECTOR 

The Beijing Electron Spectrometer, BES, is a conventional cylindrical magnetic detector 
that is coaxial with the BEPC colliding e~^e~ beams. It is described in detail in Ref. [|5|. A 
four-layer central drift chamber (CDC) surrounding the beampipe provides trigger informa- 
tion. Outside the CDC, the forty-layer main drift chamber (MDC) provides tracking and 
energy-loss (dE/dx) information on charged tracks over 85% of the total solid angle. The 
momentum resolution is (Tp/p = 0.017a/1 ± p'^ {p in GeV/c), and the dE/dx resolution for 
hadron tracks for this data sample is ~ 9%. An array of 48 scintillation counters surround- 
ing the MDC provides measurements of the time-of-flight (TOP) of charged tracks with a 
resolution of ~ 450 ps for hadrons. Outside the TOP system, a 12 radiation length lead-gas 
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barrel shower counter (BSC), operating in self- quenching streamer mode, measures the en- 
ergies of electrons and photons over 80% of the total solid angle. The energy resolution is 
E = 0.22/V^ {E in GeV). Surrounding the BSC is a solenoidal magnet that provides a 
0.4 Tesla magnetic field in the central tracking region of the detector. Three double layers of 
proportional chambers instrument the magnet flux return (MUID) and are used to identify 
muons of momentum greater than 0.5 GeV/c. Endcap time-of-flight and shower counters 
extend coverage to the forward and backward regions. 



III. TECHNIQUE 

Our measurement is based on a data sample corresponding to an integrated luminosity of 
about 6.1 pb~^ accumulated at the ^ji^S) resonance. The ^ji^S) is a copious source of J/ip 
decays: the branching fraction ip{2S) -k^-k^ J /ip = 0.324 ± 0.026 |^, is the largest single 
ipi^S) decay channel. We determine the J ftp leptonic branching fraction from a comparison 
of the exclusive and inclusive processes: 

and anything (//) 



The J/ip leptonic branching fraction is determined from the relation 



B{J/ij 1+1-) 



where Nf"^ and Nj^^ are observed numbers of events for processes I and II, and and Sj/^ 
are the respective acceptances. 



IV. EVENT SELECTION 

For both processes I and II, we use only runs of good quality. We require at least one 
pair of oppositely charged candidate pion tracks that each satisfy the following criteria: 

1. < 0.5 GeV/c, where is the pion momentum. 

2. P-nxy > 0.1 GeV/c, where P^xy is the momentum of the pion transverse to the beam 
direction. This removes tracks that circle in the Main Drift Chamber. 

3. I cos^ttI < 0.75. Here 6.,^ is the polar angle of the vr in the laboratory system. 

4. cos^^TTTT < 0.9. 9t,t, is the laboratory angle between the 7r+ and vr". This cut is used to 
eliminate contamination from misidentified e~^e~ pairs from 7 conversions, as shown 
in Fig. |I]. 

The invariant mass recoiling against the candidate tt+tt" pair, rrf^+'^i = [{m^^2S) ~ Ej^+ — 
E^-y — (p^+ + p,r-)^]^! is required to be in the range 3.0 < Tif^+'^i < 3.2 GeV/c^. 
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FIG. 1. The distribution of the cosine of the angle between the 7r+ and vr", cosOj^t^, in 
ip' — > vr+vr" J/V', J/ip ^ e^e~ events (dots with error bars) compared with the distribution of 
Monte Carlo data (histogram). The data has background near cos0,r7r = 1. 

A. J/jP £+£- 

For leptonic decay candidate events (process-I), the number of charged tracks is required 
to be at least four with a 4-track combination of net charge zero [0 . Lepton pair candidates 
must satisfy the following selection criteria: 

1. Pi > 0.5 GeV/c. Here Pi is the three-momenta of the candidate lepton track. 

2. Pi+ > 1.3 GeV/c or Pi- > 1.3 GeV/c or {Pi+ + Pi-) > 2.4 GeV/c. This cut selects 
events consistent with J/ip decay, while rejecting background. 

3. |cos6'e| < 0.75, |cos6'^| < 0.60. Here 9e and are the polar angles of the electron 
and muon, respectively. This cut ensures that electrons are contained in the BSC and 
muons in the MUID system. 

4. cos^^;'^™- < —0.975, where ^^^TJ- is the angle between the two leptons in the rest frame 
of the J/ip. 

5. For e~^e~ candidate pairs: SCE^ and SCE_ > 0.6 GeV/c, where SCE is the energy 
deposited in the BSC, or, if one of the tracks goes through a BSC rib or has Pi < 0.8 
GeV/c, the dE/dx information of both tracks in the MDC must be consistent with 
that expected for electrons. The rib region of the BSC is not used because the Monte 
Carlo does not model the energy deposition well in this region. 

6. For pair candidates at least one track must have N^"^^ > 1, where A^'^** is the 
number of MUID layers with matched hits and ranges from to 3. If only one track is 
identified in this fashion, then the invariant mass of the /i/i pair must also be within 
250 MeV/c^ of the J/ip mass. 

Fig. 2a shows the mJI^+'^L distribution for the ip' — > ti^-k'J/iI)^ J/ip events. 
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FIG. 2. a) Number of events versus m^*!f°*i, the mass recoiling against the two vr's, for 
ip{2S) TT^TT^ J /ip, J 1^ ^ events. The histogram is data, and the smooth curve is a BWG 
("signal") plus a third order polynomial ("background"), b) Number of events versus rrf^+'^L for 
inclusive events. The histogram is data, and the smooth curve is a BWG (signal) with parameters 
determined from a) plus a fourth order polynomial (background). 

B. Fitting 

The number of process-II events {J /ip — * anything) is determined from a fit to the tt+tt" 
recoil mass spectrum, using a. J /ip hne shape that is determined from the recoil mass spec- 
trum for the leptonic decays (Fig. 2a). This is fit with a Breit-Wigner function folded with 
a Gaussian (BWG) and a third order background polynomial [Q. The parameters thus ob- 
tained for the BWG fit are then used together with a fourth order polynomial background 
to fit the 7r+7r~ recoil mass for the inclusive J /ip decays, as shown in Fig. 2b. The resulting 
number of inclusive decays is Nf"^ = 530423 ± 1270, where the error is the statistical uncer- 
tainty combined with the uncertainty in the fitting procedure. The same BWG parameters 
and a third order background are used to obtain the number of e+e" and events sep- 

arately and yield a total of 18118 ± 150 J/tp —>■ e^e~ events and 14611 ± 134 J/ip A*'''/^" 
events. The results are summarized in Table |. 

TABLE I. Summary of Results. 





J /ip anything 


J/V' e+e 






530423 ± 1270 


18118 ± 150 


14611 ± 134 




44.67% 


25.85% 


21.07% 


j^qobs I^MC 


1.1876 X 10*^ 


70089 


69345 
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V. ACCEPTANCE 



The acceptances are obtained from Monte Carlo simulations. According to ref. the 
orbital angular momenta between the tt+tt" system and the J/V', as well as that between 
the TT"*" and 7i~ is zero, and the tt+tt^ mass, mx, distribution is 



da 
dmx 



oc (Phase Space) x {M^ - 4m^ 



2\2 



OC {Ml - Amlf{Ml - Amlf'^ x 



,2 \2 



4m 



2 11/2 



We generate process-I as sequential two-body decays ip{2S) X + J /ip, X tt+tt^, and 
J/il) ^ . Isotropic angular distributions are used for the J /if) in the laboratory and for 
the charged pions in the X rest frame |10[ . Leptons are generated with a 1 + cos^ 9 angular 
distribution in the J/ip rest frame and with order final state radiative corrections [jlT . 
The TT"*" and vr" decay in the detector according to the PDG |^ life time and branching 
ratios. Initial state radiation is not included in the ip{2S) generation for data taken at the 
ip{2S) peak energy. See Fig. ^ for a comparison of the mx distribution for data and Monte 
Carlo generated data. 




0.3 0.35 0.4 0.45 0.5 0.55 0.6 

mx (GeV/c^) 

FIG. 3. Number of events versus mx- The points with error bars are data, and the histogram 
is Monte Carlo data. 

Monte Carlo samples of about eight times the number of events produced in our data 
samples via process-I are generated. After application of the same selection criteria and 
fitting procedure as used for the data, we get Se = 25.85% and = 21.07%. 

For process-II, the nn acceptance, ej/^, depends on the charged particle multiplicity 
produced in the J/ ip decay. The acceptances obtained from generating different multiplicity 
J/ ip decay events are listed in Table y. Also shown is the multiplicity distribution obtained 
from our data [0. Using this distribution and the acceptances, we obtain ej/^ = (44.67 ± 
0.20)%, where the error comes from varying the multiplicity values over the range of values 
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reported by other experiments |T3|, as well as the Monte Carlo statistical uncertainties. The 
acceptance ejj^ is not sensitive to the number of vr^'s or K°'s accompanying the charged vr's 
in the J/V' decay channel or to replacing the charged vr's with charged kaons. 



TABLE II. Acceptance (vrvr) for different charged multiplicities. 



Charged 


MC decay 


Percent 


Acceptance (%) 


Multiplicity 


channel 


of J/V' 







all neutrals 


0.5 


48.34 


2 




35.6 


46.47 


4 


2(7r+7r-) 


40.9 


44.26 


6 


3(7r+7r-) 


17.7 


43.03 


8 


4(7r+7r~) 


5.3 


41.32 



VI. ERROR ESTIMATION 

jobs 



The statistical branching fraction errors, propagated from the statistical errors on A^";"' 
and A^f ^ are an = 0.0005 and (Xn = 0.0006. 



The effect of changing cuts has been studied. The results are shown in Table |T]. 
Other systematic contributions come from the acceptance uncertainties, particle identifi- 
cation (PID) uncertainties, and the fitting method. For ee events, electron particle identifi- 
cation backgrounds and efficiencies can be measured by comparing the BSC results and the 
dE/dx information from the MDC for events outside the BSC ribs. The backgrounds found 
and the difference in efficiency between that determined from the data and that predicted 
by the Monte Carlo are small and are used in the estimate of the electron PID systematic er- 
ror. The PID systematic error for the /x/x events is obtained by determining the background 
allowed by the invariant mass cut for events where only one track is identified by 
the MUID system and the estimated efficiency of this cut. The systematic error associated 
with the fitting procedure is determined by using an alternative fitting method. The total 
systematic error is taken as the sum in quadrature of all the individual systematic errors. 
The relative systematic errors on both the J/ip ^ e^e~ and ii^ ii~ branching fractions are 
1.7%. 

Contamination from hadronic events has been checked by using a sample of kinematically 
selected ip{2S) ti^ti^ J /ip, J /ip —>■ pvr events. The number of these events satisfying 
the cuts used in this analysis was negligible. In addition, a simulation of the potential 
background process ip{2S) —>■ rjj/ip with rj going to 7r+7r~7r° or 7r+7r~7 indicates a negligible 
contribution. 
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TABLE III. Branching Ratio Systematic Errors (%) 



Variable 


Variation 


i3(ee) 


Biuu) 

^ KH'H'J 


ros fl rnt 


75 ^ 70 


0.21 


0.74 


COS B cut 


9 ^ 85 


0.07 


0.07 


P cut 


tiTTTi nfr 


0.46 


0.56 






05 




1 COS Qe 1 cut 


0.75 0.7 


0.19 


0.10 


COS Ql^i- cuts 


turn off 


0.46 


0.28 


SCE cut 


0.6 ^ 0.7 GeV 


0.18 




or cut 


> 1.3 ^> 1.4 GeV/c 


0.30 


0.01 


and P;^ cut 


> 0.5 ^> 0.8 GeV/c 


0.89 


0.09 


Use only best tracks 




0.14 


0.43 


Fitting Method 




0.85 


1.21 


Efficiency Uncertainty 




0.41 


0.45 


PID Uncertainties 




0.84 


0.17 


Combined 




1.69 


1.71 



VII. RESULTS AND DISCUSSION 



The numbers of events obtained are summarized in Table | fT^, and the final results for 
the branching fractions are: 

B{J/tp e+e") = (5.90 ± 0.05 ± 0.10)% 

and 

B{J/ip fi+fi-) = (5.84 ± 0.06 ± 0.10)%. 
The close equality of and B^ is a verification of e-fi universality: 

B{J/iIj e+e-) 



1.011 ±0.013 ±0.016. 



Assuming B^ = B^ we find a combined leptonic branching fraction of 

B{J/i) tr) = (5.87 ± 0.04 ± 0.09)% 

Our results are compared with previous experiments in Table They are consistent 
with and improve on the precision of the MARK-III measurement |^ . They are also consis- 
tent with BES results determined from e~^e~ cross section measurements in the vicinity of 
J/i/j resonance If we combine the values of B{ll), we obtain a new world average 

B{ll) = (5.894 ±0 

which has an error about half that in the 1996 PDG 
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TABLE IV. Comparison with other experiments 



Experiment B{ee){%) B{^iii){%) 

Mark I § 6.9 ± 0.9 6.9 ± 0.9 

Mark III g 5.92 ± 0.15 ± 0.20 5.90 ± 0.15 ± 0.19 

BES Scan [§ 6.09 ± 0.33 6.08 ± 0.33 

THIS EXPERIMENT 5.90 ± 0.05 ± 0.10 5.84 ± 0.06 ± 0.10 



The experimental ratio of the quarkonium annihilation rates 

T [quarkonium ggg) 



r {quarkonium fi^fi 



allows the determination of the strong coupling constant or, equivalently, the QCD scale 
parameter A^^^ p|JT6[]. Following the notation of ref. |T^, one obtains 

(1 + Ri)R2 = [1 - 2B{J/iIj l+r) 

- B{J/i; 7* ^ qq) 

- B{J/^ ^ 7r?,)]/i?(J/^ ^ /+/-) 



where, 



^ r(J/^ igg) ^ i6a -30— 
^ r(J/V' ggg) 5a, ' tt 

T{J/ij ^ 5(7r2 - 9)« 
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and 



r(cc — > qqq) 
7^= J. - =0-31 ±0.03. 



The factor 7^ is the reduction factor of the three gluon decay caused by the finite size effect 
in the matrix element. For we use the theoretically calculated rather than the measured 



value because of the large error associated with the latter [17 



Combining the above J/ip decay parameters and the values for B{J/^jj — > 7* — ;> qq) and 
B{J/ip jric) listed in PDG [H, we obtain 

as{mc = l.5GeV/c^) = 0.28 ± 0.01 and 

Ag^ = (209 ± 21)MeV, 

where only the experimental errors are included. Our results are in good agreement with 
the PDG § values. 

We would like to thank the staff of BEPC accelerator and the IHEP Computing Center 
for their efforts. 
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